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nonlinearly on the fermion fields. In this framework the fermion mass 
term does not break chiral symmetry This property allows us to use 
the Wilson term to remove the doubler fermions while maintaining ex- 
act chiral symmetry on the lattice. Our lattice formulation enables us to 
address non-perturbative questions in effective field theories of baryons 
interacting with pions and in models involving constituent quarks in- 
teracting with pions and gluons. We show that a system containing a 
non-zero density of static baryons interacting with pions can be studied 
on the lattice without encountering complex action problems. In our 
formulation one can also decide non-perturbatively if the chiral quark 
model of Georgi and Manohar provides an appropriate low-energy de- 
scription of QCD. If so, one could understand why the non-relativistic 
quark model works. 
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1 Introduction 



Realizing chiral symmetry on the lattice is a problem with a long history. In 1975 
Wilson decided to break chiral symmetry explicitly in order to remove the unwanted 
doubler fermions [1]. As a result, recovering chiral symmetry with Wilson fermions 
requires fine-tuning as well as taking the continuum limit. In 1981 it was proved by 
Nielsen and Ninomiya [2] and in a different way by Friedan [3] that removing the 
doubler fermions while maintaining locality is possible only if one explicitly breaks 
the chiral symmetry of the continuum theory. In 1992, using a Wilson term in 4+1 
dimensions, Kaplan constructed domain wall fermions in order to realize chiral sym- 
metry on the lattice [4]. At the end of the same year, Narayanan and Neuberger 
introduced the related overlap formulation of chiral lattice fermions [5]. Already 
in 1982 Ginsparg and Wilson [6] had discussed how to maintain a remnant of chi- 
ral symmetry on the lattice. This paper was rediscovered by Hasenfratz in 1997 
in relation to the perfect action approach to QCD [7]. In particular, perfect lattice 
fermions obey the by now well-known Ginsparg- Wilson relation. It was immediately 
realized by Neuberger that the same is true for overlap fermions [8] . In 1998 Liischer 
showed that lattice fermion actions that obey the Ginsparg- Wilson relation realize 
a new version of chiral symmetry that is natural for theories on the lattice and co- 
incides with the usual one in the continuum limit [9]. This development culminated 
in the non-perturbative construction of chiral lattice gauge theories [10]. 

In this paper we consider another aspect of lattice chiral symmetry, namely how 
it can be realized nonlinearly. Nonlinear realizations of symmetries are well under- 
stood in the continuum in the context of low-energy effective theories describing 
spontaneous symmetry breaking. Goldstone's theorem tells us that the spontaneous 
breakdown of a continuous global symmetry G to an unbroken subgroup H implies 
the presence of massless particles: the number of Goldstone boson fields is given by 
the difference between the number of generators of G and H [11]. The chiral sym- 
metry group G = SU(N f ) L <g> SU(N f ) R ®U(1) B of QCD with N f massless quark 
flavors is spontaneously broken to the vector subgroup H = SU(Nf) L=R <g> U{1) B . 
Hence there are Nj — 1 massless Goldstone bosons: the three pions for Nf = 2 and, 
in addition, the four kaons, and the 77- meson for Nf = 3. The theoretical framework 
for describing Goldstone bosons with effective fields was developed by Weinberg [12] 
and by Coleman, Wess, and Zumino [13], partly in collaboration with Callan [14]. 
As a result, Goldstone bosons are described by fields in the coset space G/H and 
the group G is nonlinearly realized on those fields. In QCD the pions, the kaons, 
and the ?7-meson are accordingly described by matrix valued fields in the coset space 
G/H = SU(N f ) L ® SU(N f ) R ® U(l) B /[SU(N f ) L=R <g> U(1) B ] = SU(N f ). Notably, 
the pion effective theory developed by Weinberg [15] has a systematic low-energy 
expansion in Gasser and Leutwyler's framework of chiral perturbation theory [16]. 

When matter fields are included in the effective theory of Goldstone bosons, they 
also transform nonlinearly under chiral symmetry [17]. In particular, the global sym- 
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metry G is hidden in a locally realized symmetry H with a corresponding "gauge" 
field that is constructed from the Goldstone boson field. It must be pointed out 
that there is no true gauge symmetry H, just the global symmetry G. Exploiting a 
nonlinearly realized chiral symmetry, Gasser, Sainio, and Svarc incorporated a single 
nucleon in chiral perturbation theory [18]. Systematic power counting schemes for 
non-relativistic heavy baryon chiral perturbation theory were developed by Jenkins 
and Manohar [19] and by Bernard, Kaiser, Kambor, and Meissner [20]. A more 
economical scheme that maintains explicit Lorentz invariance was introduced by 
Becher and Leutwyler [21]. 

One may ask why one would want to put low-energy effective theories on the lat- 
tice. Indeed, many questions that arise in the framework of low-energy Lagrangians 
can be solved analytically using chiral perturbation theory. Usually, such calcula- 
tions are performed in the continuum using dimensional regularization. It is certainly 
more tedious to work with the lattice regularization in chiral perturbation theory 
calculations. Still, in order to address conceptual problems, for example, related 
to power divergences, the lattice has already been used in magnon [22], pion [23], 
as well as heavy baryon chiral perturbation theory [24]. In the context of lattice 
QCD, effective field theories provide powerful tools to keep systematic errors un- 
der control [25] . Moreover, chiral Lagrangians on the lattice have been investigated 
in order to facilitate the lattice determination of the Gasser-Leutwyler coefficients 
[26, 27]. Although interesting, this is not the main motivation for the present work. 
We formulate lattice theories with nonlinearly realized chiral symmetry in order to 
address non-perturbative questions that arise in the context of low-energy effective 
theories. It is perhaps not obvious that such problems even exist. An example of 
a non-perturbative problem in the context of chiral Lagrangians concerns the rotor 
spectrum of massless QCD in a finite box. However, this problem reduces to one in 
quantum mechanics and has been solved analytically by Leutwyler [28] without any 
need for a lattice. Few-nucleon systems have been treated with effective field theory 
methods in nuclear physics. For example, the two-nucleon system was studied by 
Kaplan, Savage, and Wise [29] and the three-nucleon system was investigated by 
Bedaque, Hammer, and van Kolck [30,31]. Again, the solution of these problems 
required non-perturbative techniques but not yet a lattice. Going to larger nuclei 
and eventually to nuclear matter is a non-trivial task that requires to solve more 
complicated non-perturbative problems in the effective theory [32] . Here the lattice 
provides a natural framework in which such questions can be addressed on a solid 
theoretical basis. 

A central observation of this paper is that the Nielsen- Ninomiya theorem is nat- 
urally avoided when chiral symmetry is nonlinearly realized on the lattice. In fact, 
constituent quark masses do not break a nonlinearly realized chiral symmetry. Simi- 
larly, the Wilson term — introduced in order to remove the doubler fermions — does 
no longer break chiral symmetry explicitly. Hence, exact chiral symmetry is natu- 
rally maintained using a simple Wilson action while the fermion doubling problem 
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is still solved. This raises the question how anomalies manifest themselves in this 
context. Remarkably, in effective theories with nonlinearly realized chiral symmetry, 
both in the continuum and on the lattice, the fermion fields do not contribute to 
anomalies. Instead anomalies must be added by hand in the form of Wess-Zumino- 
Witten [33,34] and Goldstone-Wilczek [35] terms. Constructing these terms on the 
lattice while maintaining their topological features is a non-trivial problem which 
we do not address in this paper. In this respect, our construction of lattice the- 
ories with nonlinearly realized chiral symmetry is still incomplete. Consequently, 
we limit ourselves to applications where anomalies play no major role. Such ap- 
plications concern, for example, models of static baryons interacting with pions at 
non-zero chemical potential and the chiral quark model of Georgi and Manohar [36]. 

Dynamical fermion simulations in a lattice theory with nonlinearly realized chiral 
symmetry may be substantially easier than with the usual linear realization. Using 
the nonlinear realization, even in the chiral limit, the fermion mass is large due to 
spontaneous chiral symmetry breaking. Hence, the simulated dynamical fermions 
should behave algorithmically like heavy Wilson fermions while still maintaining 
manifest chiral symmetry. Using a linear realization of chiral symmetry on the 
quark fields, explicit auxiliary pion fields have already been explored in order to 
simplify QCD simulations [37,38]. It has to be mentioned that, in general, lattice 
theories with nonlinearly realized chiral symmetry have a potential sign problem 
similar to the one that arises for a single flavor of Wilson fermions. Fortunately, this 
sign problem is a lattice artifact which should disappear in the continuum limit, in 
particular, for sufficiently heavy fermions. Recently, without directly addressing the 
potential sign problem, dynamical fermion algorithms have been developed for an 
odd number of flavors of Wilson fermions [39-42]. We expect that similar methods 
are applicable to lattice fermions with nonlinearly realized chiral symmetry. Still, 
it should be pointed out that the feasibility of dynamical fermion simulations de- 
pends crucially on the severity of the sign problem, which should be investigated 
in numerical simulations. We also need to mention that additional complex ac- 
tion problems arise if pseudo- vector fermion couplings or Wess-Zumino-Witten and 
Goldstone-Wilczek terms are included. In this paper we concentrate on the theo- 
retical formulation of lattice theories with a nonlinearly realized chiral symmetry, 
while results of numerical simulations will be presented elsewhere. 

Interestingly, the lattice construction presented here enables us to simulate mass- 
less pions interacting with static baryons even at non-zero baryon chemical potential 
/i and finite temperature T, both for two and for three flavors. In contrast to full 
QCD, in the static baryon model this is possible without encountering a complex 
action problem. Although questions involving high temperatures or high baryon 
densities are beyond the applicability of a systematic low-energy expansion, this 
model can still be used to answer questions concerning the universal features of the 
critical endpoint of the chiral phase transition line in the (/x, T)-plane. Although 
one should not expect to obtain precise information on non-universal features of 
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QCD, such as the location of the critical endpoint, one may learn at least semi- 
quant it at ively how the location of the endpoint changes as a function of the strange 
quark mass. In addition, one will gain experience with numerical methods for locat- 
ing the critical endpoint in a model that has the same global symmetry properties 
as full QCD but is much simpler to simulate. 

A nonlinear realization of chiral symmetry was also used by Georgi and Manohar 
in the construction of their chiral quark model [36]. This model contains massive 
constituent quarks interacting with pions and gluons. While gluons are flavor sin- 
glets, both constituent quarks and pions transform non-trivially under the nonlin- 
early realized chiral symmetry. If this model provides an appropriate low-energy 
description for QCD, one could understand why the non-relativistic quark model 
works. However, the confinement of constituent quarks is a non-perturbative issue 
that cannot be addressed quantitatively in the continuum formulation of the model. 
In this paper we present a lattice construction of the Georgi-Manohar model with 
nonlinearly realized chiral symmetry. This non-perturbative framework will allow 
one to decide quantitatively if the Georgi-Manohar model can explain the success 
of the non-relativistic quark model. 

The rest of this paper is organized as follows. In section 2 the nonlinear real- 
ization of chiral symmetry is reviewed in the continuum, while section 3 contains 
the corresponding construction on the lattice. In both cases the relevant symme- 
tries are discussed in detail. In particular, effective lattice theories for pions and 
baryons are considered and the issue of fermion doubling is discussed. In section 4 
we discuss how low-energy effective theories can be treated non-perturbatively on 
the lattice and how continuum results can be extracted. Section 5 explains how the- 
ories with static baryons can be simulated at non-zero chemical potential without 
encountering a complex action problem. In section 6 we discuss the lattice version 
of the Georgi-Manohar model and propose numerical simulations in order to decide 
if it indeed provides a satisfactory explanation for the success of the non-relativistic 
quark model. Section 7 discusses the possible form of the phase diagram of the lat- 
tice chiral quark model and the relation of this model to Wilson's standard lattice 
QCD. Finally, section 8 contains our conclusions. 



2 Continuum Formulation of Nonlinearly 
Realized Chiral Symmetry 

In this section we review how theories with nonlinearly realized chiral symmetry are 
formulated in the continuum. The fields used in this construction are introduced 
with emphasis on their behavior under chiral symmetry transformations, charge 
conjugation, and parity. Simple examples of low-energy effective theories with non- 
linearly realized chiral symmetry are given and manifestations of anomalies are being 
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discussed. 



Consider the chiral symmetry group G = SU(Nf)L <8> SU(Nf) R ® U(1)b of 
QCD involving Nf quark flavors which breaks spontaneously to the subgroup H = 
SU(Nf) L=R <S> U(1)b- Following [13, 14] the Goldstone boson field is then given by 
U(x) = exp[2m a (V)TVF^] G G/H = SU(N f ). Here Tr (T a T b ) = 25 ab and F w is the 
pion decay constant. Under global chiral rotations L <8> R G SU(Nf) L <g> SU(Nf) R , 
this field transforms as 

U(x)' = LU(x)R t . (2.1) 
From the Goldstone boson field one derives a field u(x) G SU(Nf) as 

u(x) = U(x) 1 ' 2 . (2.2) 

In order to avoid ambiguities in taking the square-root, we proceed as follows. First, 
the field U(x) is diagonalized by a unitary transformation W(x), i.e. 

U(x) = W(x)D(x)W(x)\ (2.3) 

where D(x) is a diagonal matrix of complex phases 

D(x) = diag{exp[i^i(a;)],exp[i^2(aj)],-,exp[i^ / (a;)]}. (2.4) 

At this point, the phases ipi(x) are defined only modulo 2n. In order to fix this 
ambiguity, we demand that 

Nf 

Y^ifiix) = 0, \<pi(x) - <pj(x)\ < 2tt. (2.5) 
i=i 

This uniquely fixes the phases of the eigenvalues and defines the matrix D(x) up to 
an irrelevant permutation of its diagonal elements. We then define 

D(x) 1/2 = diag{exp[i(p 1 (x)/2],exp[i(p 2 (x)/2],...,exp[i(p Nj (x)/2]}, (2.6) 

and we obtain 

u(x) = W{x)D{x) 1/2 W{x) ] . (2.7) 

When the eigenvalues of U(x) are non-degenerate, the matrix u(x) is uniquely de- 
fined in this way and is located in the middle of the shortest geodesic connecting 
U(x) with the unit-matrix 1 in the group manifold of SU(Nf). On the other hand, 
if two eigenvalues are degenerate, the corresponding eigenvectors are not uniquely 
defined. If, in addition, the corresponding phases ipi and ifij differ by 2n, this ambi- 
guity gives rise to a coordinate singularity in the definition of u(x). In the Nf = 2 
case, the singularities arise for U(x) = —1, which corresponds to a single point in 
the SU(2) group manifold. For general Nf, the submanifold of SU(Nf) matrices 
with two degenerate eigenvalues (and with \cpi — ipj\ = 2-k) has dimension Nj — 4. 
Hence, the codimension of the singular submanifold is Nj — 1 — (Nj — 4) = 3 for all 
Nf. Consequently, in a 4-dimensional space-time the singularities in u(x) are located 
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on 1- dimensional lines. If one chooses the vacuum as U(x) = 1, the singularities 
represent the world-lines of Skyrmion centers. However, it should be noted that the 
location of the singularities varies under those chiral rotations that do not leave the 
vacuum invariant. Still, the effective actions that we will construct with u(x) are 
completely chirally invariant, even in the presence of Skyrmions or other large-field 
configurations. 

Under chiral rotations the field u(x) transforms as 

u{x)' = Lu(x)V( X y = V{x)u{x)R\ (2.8) 

The matrix V(x) depends on L and R as well as on the field U (x) and can be written 

as 

V(x) = R[R^LU( X )] 1/2 [U( X ) 1/2 )^ = LpRU(x)^ 2 U(x)^ 2 . (2.9) 

For transformations in the unbroken vector subgroup SU(Nf)L = R, L = R, the field 
V(x) reduces to the global flavor transformation V(x) = L = R. The local transfor- 
mation V(x) is a nonlinear representation of chiral symmetry. For two subsequent 
chiral transformations with L 2 L\ = L and R2R1 = R, one has 

V^x) = R 1 [R\L 1 U(x)] 1/2 Mx) 1/2 }\ V 2 (x) = Ri[R\L 2 U{x)'Yl 2 [U{x)' l/2 \\ (2.10) 

with U(x)' = L 1 U(x)R\, so that indeed 

V 2 (x)V 1 (x) = R 2 [R\L 2 L l U{x)R\Y/ 2 R l [U{xf/ 2 ^ 
= R^iRlRlL^Uix^iUix) 1 ^ 
= R[R)LU{x)\ 1 / 2 [U{x) 1 l 2 Y = V(x). (2.11) 

Now consider a Dirac spinor field ty(x) and ty(x) that transforms under the 
fundamental representation of the flavor group. This field may represent a nucleon 
(for Nf = 2) or a constituent quark (for any value of Nf). Global chiral rotations 
L <S> R € SU(Nf)i <E> SU(Nf)n can be realized nonlinearly on this field using the 
t r ansf or mat ions 

y(x)' = V(x)*(x), = ^(x)V(x)l (2.12) 

This has the form of an SU(Nf) gauge transformation despite the fact that it rep- 
resents just a global SU(Nf) L £g> SU(Nf) R symmetry. 

In order to construct a chirally invariant (i.e. SU(Nf) "gauge" invariant) La- 
grangian, one needs an SU(Nf) flavor "gauge" field. Indeed, the anti-Hermitean 
composite field 

v i*( x ) = - \u{x) ] d^u(x) + u(x)d^u(x) ] ] (2.13) 
transforms as a "gauge" field 

v^x)' = ^{V(x)u(xyL' ( d fl [Lu(x)V(xy}+V(x)u(x)R' ( d l ,[Ru(xyV(xy}} 

= V(x)[v^x) + d^V(x)l (2.14) 



8 



A Hermitean composite field is given by 

a ni x ) = l -[u(xyd^u(x) - u(x)d fl u(xy], (2.15) 

which transforms as 

a^x)' = ^{V(x)u(xytfd^{Lu(x)V(xy] - V {x)u{x)F$ d^[Ru{x) ] V {x^}} 

= V(x)a^(x)V(x)l (2.16) 

It is interesting to investigate the transformation properties of the fields *&(x), 
^(x), U(x), u(x), Vfj,(x), and a^(x) under charge conjugation and parity. Under 
charge conjugation, the fermion field transforms in the usual way 

c ^( x ) = CV{x) T , c V{x) = -<S>(x) T C-\ (2.17) 

with C~ 1 7^C = — 7 J. The Goldstone boson field U(x) transforms as 

c U(x) = U(x) T . (2.18) 

Consequently, the field u(x) also follows the same transformation law 

c u {x) = [ C U(x)\ 1 / 2 = [Uixf] 1 / 2 = [Uixf'Y = u{x) T . (2.19) 

It is easy to show that the composite fields v^(x) and a^x) transform as 

C v,(x) = v^x)*, %(x) = a^xf. (2.20) 

Under parity, the fermion field behaves again in the usual way, i.e. 

p ^{x , t) = 74* (-£, t), P ^(x, t) = (-£ , t) 74 , (2.21) 

where 4 denotes the Euclidean time direction. The Goldstone boson field U(x) 
transforms as 

p U(x,t) = U(-x,t) ] (2.22) 

so that 

p u(x, t) = [ p U(x, t)} 1 ' 2 = [U(-x, tY] 1 / 2 = [U(-x, tf 2 ^ = u(-x, t)l (2.23) 
Similarly, one finds 

p Vi(x,t) = -Vi(-x,t), p v 4 {x,t) =v 4 {-x,t), 

p aj(x, t) = ai(—x,t), p a 4: (x,t) = — a 4 (— x, t), (2.24) 

where % G {1, 2, 3} denotes a spatial direction. Hence, is a vector and a^{x) is 

an axial vector. 

9 



At this point we are prepared to write down the leading terms in the Euclidean 
action of a low-energy effective theory for nucleons (with Nf = 2) 

S\V, #,[/] = J d 4 x {AW + *7m(^ + + WA^l^a^ 

+ ^Tr[d^d^U]-^^Tr[MU* + M*U}}. (2.25) 

Here M is the nucleon mass generated by spontaneous chiral symmetry breaking, 
qa is the coupling to the isovector axial current, and (iptp) and M. = diag(m u , m^) 
are respectively the chiral condensate per flavor and the mass matrix of the current 
quarks. When ty(x) represents a constituent quark, M would be the constituent 
quark mass. In that case, gluons must also be included in order to confine the 
constituent quarks. This is what happens in the chiral quark model of Georgi and 
Manohar [36] (see section 6). 

It is remarkable that — due to the nonlinear realization of chiral symmetry — the 
fermion mass term is chirally invariant. This makes sense because the mass M arises 
from the spontaneous breakdown of chiral symmetry even in the chiral limit, i.e. 
when the current quarks are massless. In fact, the term that contains the matrix M. 
of current quark masses is the only source of explicit chiral symmetry breaking in the 
leading terms of the effective theory. It is interesting to note that the composite field 
v^x) is coupled to the fermions as a flavor "gauge" field that enters the covariant 
derivative = + v^. Like i> M , the axial vector field a M also transforms in the 
adjoint representation of the flavor group. However, it is not a "gauge" field and 
it couples to the fermions with strength qa- It must be noted that there is, of 
course, no true SU(Nf) gauge symmetry, just an SU(Nf) L <g> SU(Nf) R <g> U(1)b 
global symmetry. In particular, the "gauge" field v^x) is not a fundamental field 
and there is, for instance, no SU(Nf) Gauss law. 

Until now we have considered a fermion field ^(x) transforming in the funda- 
mental representation of SU(Nf). In QCD with three massless flavors (Nf = 3) the 
baryons transform as a flavor octet. The corresponding Dirac spinor field, denoted 
by x( x ) an d x{ x )i 1S traceless and transforms in the adjoint representation 

X (x)' = V(x) X (x)V(x)l (2.26) 

The covariant derivative is then given by 

D»x( x ) = d M x ) + M^)> X(x)}, (2.27) 

which transforms as D^x( x )' — V (x) D ^x( x )V ( X Y under the nonlinearly realized 
chiral symmetry. In this case, the Euclidean action of the effective theory for pions, 
kaons, 77-mesons, and baryons takes the form 

S[x, X,U] = J d 4 x {Tt[Mxx + Xl^pX + ^X7^75 W, x} + ^XT/^sfS' x\\ 

+ ^Tr[d^d^U]-^Tr[MU^ + M j U]}, (2.28) 
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where the matrix M. = diag(m n , ma, m s ) now also includes the strange quark mass. 
The low-energy parameters D and F are analogous to the coupling g& of the two 
flavor case. 

It is interesting to discuss how anomalies arise in effective theories with nonlin- 
early realized chiral symmetry. Let us begin with a discussion of the axial U{1)a 
symmetry of QCD. When the number of colors N c is finite, this symmetry of the 
classical action is anomalously broken by quantum effects due to instantons [43]. 
As a result, the rf- meson is not a Goldstone boson. Only in the large N c limit 
the T/'-meson becomes a pseudo-Goldstone boson [44,45]. For N c = 3, the U(1)a 
symmetry is so strongly broken that it cannot even be considered an approximate 
symmetry of QCD. In particular, the 7/-meson is so heavy that it is not included 
explicitly in the low-energy effective theory. Indeed, the effective theories that we 
have written down do not have a U(1)a symmetry at all. This follows because the 
Goldstone boson field lives in SU(Nf) and not in U(Nf). 

Let us now move on to those anomalies that arise when one gauges the elec- 
troweak interactions. In effective theories with nonlinearly realized chiral symme- 
tries anomalies manifest themselves [17] in the form of Wess-Zumino-Witten [33,34] 
and Goldstone- Wilczek terms [35]. For example, in the Nf = 2 case, Witten's 
global SU(2)l anomaly [46] manifests itself in a topological factor Signff/]^ G 
II4 [iS'Z7(2)] = {1, —1} that must be included in the path integral [47]. This ensures 
that, for odd N c , the Skyrme soliton [48] is quantized as a fermion representing the 
physical nucleon in the pure pion theory. The SU(2) l gauge invariant version of the 
Skyrme baryon number current is the Goldstone- Wilczek current [35]. In the pres- 
ence of electroweak instantons, this current is no longer conserved and reproduces 
the 't Hooft anomaly in the baryon number [43]. 

In the Nf > 3 case, the Wess-Zumino-Witten term replaces the topological factor 
Sign[C/] in the Goldstone boson theory. This is because then U^SUlNf)} = {0} and 
instead U 5 [SU(Nf)] = Z. The Wess-Zumino-Witten action has the form 

(2.29) 

where H 5 is a 5-dimensional hemi-sphere whose boundary dH 5 = S 4 is (compacti- 
fied) space-time. The physical Goldstone boson field U(x) G SU (Nf) is extended to 
a field U(x, £5) that depends on the unphysical fifth coordinate X5 G [0, 1] such that 
U(x,0) = 11 and U(x,l) = U(x). Of course, the 4-dimensional physics should be 
independent of how the field U{x,x$) is deformed into the bulk of the fifth dimen- 
sion, and should only depend on U(x). This is possible because the interpolation 
ambiguity of SVzty^] is given by 2n times a winding number in H 5 [SU(Nf)] = Z. 
Hence, when the prefactor of the Wess-Zumino-Witten term is quantized in integer 
units, the 4-dimensional physics is independent of how one interpolates the Gold- 
stone boson field into the unphysical fifth dimension. Remarkably, in the low-energy 
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theory for QCD, the quantized prefactor is the number of colors N c so that 

Z = J VUexp(-S[U})exp(iN c S WZ w[U}). (2.30) 

For Nf = 2, one can show that exp(iN c Swzw[U}) = Siga[U] Nc . The Wess-Zumino- 
Witten term contains the vertex that describes the anomalous decay of the neutral 
pion into two photons. For general N c , this process also receives contributions from 
a Goldstone-Wilczek term [49]. 

In contrast to the Nf = 2 case, for Nf > 3, the Wess-Zumino-Witten term also 
breaks an unwanted intrinsic parity symmetry P that is present in the Goldstone 
boson action S[U] but not in QCD [34]. The full parity P acts on the pseudo-scalar 
Goldstone bosons ir a (x,t) by spatial inversion accompanied by a sign-change, i.e. 
p ir a (x,t) = —7T a (—x,t), such that 

p U(x,t) = U(-x,tf. (2.31) 

The intrinsic parity P , on the other hand, leaves out the spatial inversion and takes 
the form 

p °U(x,t) = U{x,t) ] . (2.32) 

If Po were a symmetry of QCD, the number of Goldstone bosons would be conserved 
modulo two. While this is indeed the case for Nf = 2 (where Po is nothing but G- 
parity [50]), for Nf > 3 intrinsic parity is not a symmetry of QCD. For example, the 
(/•-meson can decay both into two kaons or into three pions. However, the Goldstone 
boson action S[U] is indeed invariant under P , 

S[ P °U] = S[U^} = S[U}. (2.33) 

Hence it has more symmetry than the underlying QCD action. The Wess-Zumino- 
Witten action, on the other hand, is odd under P , 

Swzw[ Po U] = S WZW [U^] = -S WZW [U], (2.34) 

and thus reduces the symmetry of the effective theory to the one of QCD. 

The next step is to ask how anomalies are represented when explicit baryon 
fields are present in the low-energy effective theory. In the presence of electroweak 
SU(2) L ® U(1) Y gauge fields W^x) = igW^(x)T a and B^x) = ig'B^T 3 one 
obtains 

v »( x ) = \{ u { x ) ] \ d n + Wp{x)\u{x) + u(x)[d„ + B^x)]u(xy}, 

= T^) 1 !^ + W^x)]u(x) - u(x)[d^ + B^x)]u( X y}, (2.35) 
which still transform as 

v^x)' = V{x)[v»{x) + d^Vixy, a^x)' = V {x)a^{x)V {x)\ (2.36) 
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Hence, the dynamical gauge fields W^x) and B^x) are incorporated in the "gauge" 
field v^{x) constructed from the pion field U(x). By varying the resulting action 
with respect to W^(x) and B^(x), one obtains the left- and right-handed fermion 
currents 

j; L (x) = nxMx)^ 1 ~ 9 2 Al5 T a u(x)V(x), 

= *(x)u(xhJ-^^T*u(xMx). (2-37) 

Since fermions with a nonlinearly realized chiral symmetry do not contribute to 
anomalies [17,51], even in the presence of fermion fields, the anomalies are still 
contained in Goldstone-Wilczek and Wess-Zumino-Witten terms. 

With explicit baryon fields being present in addition to the Skyrme solitons in the 
pion field, one may wonder if baryons are doubly counted in this approach. First of 
all, in the power counting scheme of chiral perturbation theory one cannot address 
questions about Skyrmions because all contributions from higher order terms are 
then equally important. In a non-perturbative approach to the effective theory, one 
can ask if, for example, in the sector with one baryon, the perturbative vacuum 
configuration of the pion field becomes unstable against Skyrmion formation. If so, 
the chiral power counting scheme is not applicable and the systematic low-energy 
approach fails. Hence, in the low-energy effective theory, one is restricted to the 
trivial topological sector and the issue of double counting does not arise. Since 
Skyrmions have no place in the systematic low-energy expansion, one cannot use 
the Goldstone-Wilczek current to describe baryon number violating processes in the 
framework of chiral perturbation theory. This should not be too surprising. Baryon 
decay — as rare as it may be — is a violent event that releases a lot of energy and 
can hence not be understood within a low-energy effective theory. 



3 Lattice Formulation of Nonlinearly Realized 
Chiral Symmetry 

Let us now construct theories with a nonlinearly realized chiral symmetry on the 
lattice. In analogy to the previous section, we set up the lattice fields needed for this 
construction such that the correct symmetry properties are respected. With these 
fields, we formulate lattice actions of low-energy effective theories in which chiral 
symmetry is nonlinearly realized. In particular, we discuss how fermion doublers 
can be removed using a Wilson term without explicitly breaking chiral symmetry 
and how the Nielsen-Ninomiya theorem is avoided. The issue of anomalies is also 
addressed in the lattice formulation. 

The Goldstone boson field U x G SU(Nf) lives on the sites x of a 4-dimensional 
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hyper-cubic lattice and it transforms under global chiral rotations as 

U' x = LU X R\ (3.1) 

As in the continuum, the field u x G SU(Nf) is constructed as u x = U^ 2 and it 
transforms as 

u' x = Lu x Vl = V x u x R ] . (3.2) 

Again, the field V x depends on the global chiral transformations L and R as well as 
on the field U x and is given by 

V x = R(R^LU x y/ 2 (Ul/y = L{l)RUl) 1 l 2 U 1 J 2 . (3.3) 

The continuum results for two subsequent chiral transformations given in eqs. (2.10) 
and (2.11) remain completely unchanged on the lattice. Also the discussion of the 
coordinate singularities still applies. However, it should be noted that the singulari- 
ties generically fall between lattice points. Still, using an appropriate interpolation, 
it is possible to detect these topological objects even on the lattice. 

Next we consider the fermion field ^f x and ^f x which again lives on the lattice 
sites x. As in the continuum, global chiral rotations are nonlinearly realized on this 
field such that 

% = V x * x , % = ^ x Vl (3.4) 

This field transforms in the fundamental representation of SU{Nf) and can hence 
represent a nucleon (for Nf = 2) or a constituent quark (for any value of Nf). 
For Nf = 3, the lattice baryon field Xx and x x a l so uves on the lattice sites but 
transforms in the adjoint representation 

X' x = V xXx Vl x' x = V x x x Vl (3.5) 



We now proceed to the construction of V Xtll G SU {Nf) — the lattice analog of 
the continuum flavor "gauge" field v^x). It should be noted that V XtjX is a flavor 
parallel transporter along a lattice link in the group SU(Nf), while the continuum 
flavor "gauge" field is in the algebra su(Nf). In analogy to the continuum expression 
(2.13) for Vn(x) : we construct the lattice object 

V Xtli = ^ + u x ul +fl ]. (3.6) 

Here fi is the unit-vector in the /z-direction. By construction, V Xjli transforms as a 
parallel transporter under the nonlinearly realized chiral symmetry, i.e. 

K, = V*V*X+n- (3-7) 

In the classical continuum limit, V x ^ = exp[t> M (x + fi/2)}. However, at finite lattice 
spacing, V Xifl is in general not an element of SU(Nf), just a complex Nf x Nf 
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matrix in the group GL(Nf). Hence, it is natural to project a group- valued parallel 

transporter V X:IX G SU(Nf) out of V x ^ by performing a GL(Nf)/SU(Nf) coset 
decomposition 

V X:IX = H x ^ exp(i<p Xili /N f )V Xili . (3.8) 

Here 

H** = (V^Vl) 1/2 ( 3 - 9 ) 
is a positive Hermitean matrix that transforms as 

H' x ^ = V x H Xtll Vl (3.10) 

and exp(i<£> Xj/1 ) (with ip x ^ e ]— 7r, 7r]) is the phase of the determinant of \4 iM . Then, 
by construction, 

V x „ = eM-^jN f )H-lV x ,„ (3.11) 
is indeed in SU(Nf) and transforms as 

= V x V x y x ^. (3.12) 

We also want to construct a lattice version of the continuum field a^x) defined 
in (2.15). For this purpose, we first consider 

Ax,i* = \ [ulux+fi - u x u x+jl } . (3.13) 

While in the continuum a^x)' = V(x)a fl (x)V(xY , the lattice field A X:fl transforms 
as 

Also, in contrast to the continuum field a^x), the lattice field A Xtfl is in general 
neither traceless nor Hermitean. To cure these problems, it is natural to introduce 
the field 

= \[A X ^ + V^AlJ - -±-Tt[A x ^ + V x ,„AtJl, (3.15) 

which, by construction, is indeed traceless and Hermitean and which transforms as 

A L x ; = V x Al i yl (3.16) 

with the matrix V x located at the site x on the left end of the link (x, /i). Similarly, 
we define the object 

= \K,K, + Al,V x J - ^Tr[Vl»A x ,„ + At^jl, (3.17) 

which transforms as 

^x,/i = ^x+p,A X4J Y x+(l , (3.18) 
15 



with the matrix V x+ ^ located at the site x + fi on the right end of the link (x, //). It 
should be noted that A£ and are not independent, but are related by parallel 
transport, i.e. 

^, = ^%- (3-19) 

Let us now consider the symmetry properties of the lattice fields under charge 
conjugation and parity. The site fields *& x , ty x , Xx, X x , U x , and u x transform exactly 
as their continuum analogs and will not be discussed again. Instead, we concentrate 
on the fields V x ^, A% , and A R . Under charge conjugation, they transform as 

C V — v* c A L A L * c A R - A R * (iVft} 
Similarly, under parity transformations, one finds 



P V- - 

Vx,t,i — 




, P Kr it) 4 - 


= V-x,t,4 j 


P A L - ■ = 


A R . ',' 


PaL 

> ^x,t,4 


- A L 


p A R t ■ = 


A L 

■™-—x,t,i ' 


PaR 

A x,t,4 — 


aR 

A -x,t,4 i 



(3.21) 

where % G {1,2,3} denotes a spatial direction and 4 again denotes the Euclidean 
time direction. Hence, as in the continuum, V Xyil is a vector and and A^ define 
an axial vector. 

We now use the fields introduced above to construct lattice actions for baryon 
effective theories. We start out with the Nf = 2 case. In close analogy to the 
continuum action (2.25), we define the lattice action as 

x ^ x,n 

x,ft 



x,n 
2 



(3.22) 

Similarly, for iVj = 3, (where the baryon field transforms in the adjoint representa- 
tion) the lattice action corresponding to (2.28) is given by 

S[x, X,U] = M E ^XxXx + \Y, MXxl^x^Xx+nVl^ - Xx+innVx^XxVx^) 

X x,\x 

iD 

+ -y E Tr fe7M75{^ lAt , Xx} + Xx+aWsi^, Xx+f,}) 

x,n 
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+ -7T ^(XxiniMx,^ Xx] + Xx+awl^, xx+fi]) 

T' 

+ « J2 Tr ( 2 X x Xx - XxVx.nXx+fiVln - Xx+pYlvXxV^) 

x,ft 

- jE nU x Ux + , + Ut + ,U X ] - M £ TV[MC/t + ^tf/j. (3.23) 

The Wilson term proportional to r removes the doubler fermions. In standard 
lattice QCD this term breaks chiral symmetry explicitly. Remarkably, when chiral 
symmetry is nonlinearly realized, not only the fermion mass term (proportional to 
M) but also the Wilson term is chirally invariant. The only source of explicit chiral 
symmetry breaking is the current quark mass matrix M.. 

It is interesting to ask how the Nielsen-Ninomiya theorem [2] has been avoided. 
Clearly, the actions (3.22) and (3.23) are local. For Nf = 2, the corresponding Dirac 
operator is given by 

D xy [U\ = M5 X)V + ^(^V x ^S x+(i! y - 7mK-a,/A-£,s/) 

+ 1 -y to^, A* + -/<-•-.•>''' ,,/<-,) 

+ ^(2Sx,y ~ V x,^x+fi,y - Vl-tJx-^y)- (3.24) 

The Nielsen-Ninomiya theorem assumes that the Dirac operator anticommutes with 
75. This is not the case when chiral symmetry is nonlinearly realized, and hence 
one of the basic assumptions of the no-go theorem is not satisfied. Interestingly, 
Ginsparg- Wilson fermions evade the Nielsen-Ninomiya theorem by violating the 
same assumption. Of course, the nonlinear realization of chiral symmetry requires 
an explicit pion field which is not present in the fundamental QCD Lagrangian. 

Let us again discuss the issue of anomalies. On the lattice, anomalies may 
receive contributions from the doubler fermions. Hence, removing the doublers is 
intimately related to anomaly matching. It is straightforward to include electroweak 
gauge fields in the lattice theory as well. For example, the parallel transporters 
W Xitl G SU(2)l of the electroweak gauge field transform as 

K* = L * W ** L U- (3-25) 
Together with the U(l)y gauge field B x ^, the field W X:fl is incorporated in the field 

Vx,p = \ HW x ^u x+fl + u x B x ^ul +ji \ , (3.26) 

which is again projected on a parallel transporter V Xjll G SU(2) by a GL(2)/ SU(2) 
coset decomposition. Also the fields and A^ are constructed as in (3.15) and 
(3.17), now using 

A x ,» = | [ulW x ^u x+fl - u x B x>lt ul +j3 \ . (3.27) 
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The lattice fermion action is constructed as in the ungauged case (3.22) and is now 
SU(2) L ® U(1) Y gauge invariant. Remarkably, not only the lattice fermion action 
but also the lattice fermion measure is gauge invariant. At first sight, this seems 
to be a severe problem because gauge invariance of both the fermion action and 
the fermion measure implies that lattice fermions with nonlinearly realized chiral 
symmetry do not contribute to anomalies. One might even suspect that the doubler 
fermions have not been properly removed and thus have canceled the anomalies of 
the physical fermions. Fortunately, this is not the case. Indeed, as we have already 
seen in the continuum theory, fermions with a nonlinearly realized chiral symmetry 
do not contribute to anomalies. Instead, even in the presence of fermion fields, the 
anomalies are still contained in Goldstone-Wilczek and Wess-Zumino-Witten terms. 
Consequently, these terms should also be included in the lattice action. It is an inter- 
esting problem to construct lattice versions of Goldstone-Wilczek and Wess-Zumino- 
Witten terms maintaining the topological properties of the continuum theory. Since 
these terms do not play a central role in the applications discussed in the remainder 
of the paper, we defer this problem to future investigations. It should be noted 
that the low-energy effective theory without the addition of Goldstone-Wilczek and 
Wess-Zumino-Witten terms, although consistent, does not correctly describe the 
low-energy physics of the underlying microscopic theory with the anomalies. We 
also would like to point out that the topological terms give rise to complex action 
problems which may severely affect the efficiency of numerical simulations. 



4 Non-Perturbative Treatment of Low-Energy 
Effective Theories on the Lattice 

In this section we describe how non-perturbative calculations can be performed in 
low-energy effective theories on the lattice. In particular, we discuss how continuum 
physics can be extracted. As we have discussed in the introduction, non-perturbative 
problems indeed exist [28-31] but have so far not required a lattice. This is likely 
to change when more complicated effective theories describing, for example, nuclear 
matter [32] are studied. 

For simplicity, we first describe how the pure pion low-energy effective theory 
can be treated on the lattice. It should be noted that, at present, there is no need 
for such calculations. The only non-perturbative problem we are aware of in this 
context is the SU(Nf) chiral rotor which has been solved analytically [28]. The 
leading term in the chiral Lagrangian for pions reads 

S[U] = I d 4 x {^fTr[d^d^U} - iM Tr [.Mf/t + M ] U}\. (4.1) 
J 4 2 
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It is straightforward to discretize this theory. A simple lattice action reads 

S[U] = -^ETr[^^ + , + Ul +fl U x ] - ^ETr[-Mf/J + M*U X ]. (4.2) 

Of course, the bare lattice parameters F w and (tpi/j) get renormalized. A basic state- 
ment of the low-energy effective theory is that, to leading order, the dynamics of 
pions depends only on the renormalized parameters F* cn and (ipip) ren . In the chiral 
limit, M. = 0, the above lattice action has a phase transition at a critical bare 
coupling F wc . For Nf = 2, the transition is second order, while for Nf > 3, it is 
expected to be first order. For F n < F nc , one is in the symmetric phase without 
Goldstone bosons, while for F n > F nc , one is in the broken phase in which the 
Goldstone bosons govern the low-energy physics. A non-perturbative calculation in 
the low-energy effective theory can be performed anywhere in the broken phase. In 
particular, there is no need to approach a second order phase transition in order to 
extract continuum physics. This is because Goldstone bosons are naturally massless 
without fine-tuning. However, the renormalized parameter F* en will in general stay 
at the cut-off scale. Still, the low-energy physics (which takes place far below the lat- 
tice cut-off) is universal and depends on the details of the lattice action only through 
F* en . In practical lattice calculations, one can choose any convenient bare coupling 
F n > F nc in the broken phase and determine the corresponding renormalized cou- 
pling F* en . Then one can choose A4(ip^p) such that the physical pion mass M^ en is 
correctly reproduced (in units of F* en ). After the low-energy parameters have been 
fixed to experiment, any non-perturbative calculation of a low-energy phenomenon 
can be performed in the lattice effective theory. Remarkably, at leading order, no 
lattice artifacts arise although F* en is at the cut-off scale. Of course, in the lattice 
theory Lorentz invariance is reduced to discrete translations and rotations. Hence, 
there are terms in the low-energy theory of the lattice model which are forbidden 
by Lorentz symmetry but not by the discrete rotation subgroup. Such terms are 
absent at leading order and arise only in the higher orders of the low-energy ex- 
pansion. Their prefactors are additional low-energy parameters (besides the usual 
Gasser-Leutwyler coefficients) which must be tuned such that, to the given order of 
the momentum expansion, the low-energy theory is free of lattice artifacts. 

Non-perturbative calculations in the baryon effective theory can be performed in 
a similar way. The effective theory correctly describes the low-energy pion dynamics 
in the sector with one baryon. In this case, the bare parameters M and qa (for 
Nf = 2) or F and D (for Nf — 3), as well as F w and A4(ipip) need to be tuned such 
that the corresponding renormalized parameters match the experiment. Just like 
F«m, now a j go baryon mass M rcn stays at the cut-off scale. This is no problem 
because at low energies the results of the effective theory are again universal and 
depend on the lattice details only through the values of the low-energy constants. 
There is no doubling problem as long as the mass of the doubler fermions is outside 
the energy range of the effective theory. In particular, there is no need to approach a 
second order phase transition in order to make the physical baryon much lighter than 
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the doubler fermions. Instead, by including higher order terms, one can push the 
doubler fermions to higher energies and thus extend the validity range of the effective 
theory. More practical aspects of simulations of lattice theories with nonlinearly 
realized chiral symmetries will be discussed in the following sections. 



5 Static Baryons at Non-Zero Chemical Potential 

Nonlinear realizations of lattice chiral symmetry may also be useful beyond the strict 
validity range of low-energy effective theories. In this section we discuss baryons in 
the static limit, M — > oo, where interesting simplifications arise. In particular, one 
can then perform simulations at non-zero baryon chemical potential without encoun- 
tering a complex action problem. We first discuss the problem in the continuum 
and then formulate the static baryon model on the lattice. Within this model we 
propose exploratory studies of the QCD phase diagram at non-zero baryon chemical 
potential. 

In the continuum formulation, static baryons at an undetermined position x 
(with = 0) are described by the spatial integral of the Polyakov loop 

$[[/] = J d 3 x Tr Vex.p[J^ dt v 4 (x,t)] (5.1) 

of the flavor "gauge" field v^{x) which is a function of u(x) = U(x) 1 ^ 2 . Here V 
denotes path ordering and (3 = 1/T is the inverse temperature. The partition 
function of pions in the background of a single static baryon then takes the form 

Z B = jvu $[[/] exp(-S[U]) exp(-PM), (5.2) 

where the baryon mass M will ultimately be sent to infinity and S[U] is the pure 
pion action given in (4.1). Similarly, the pion partition function in the presence of 
a single static antibaryon at an undetermined position is given by 

% = JvU $[[/]* exp(-S[U]) exp(-/3M). (5.3) 

In the next step, we consider a system of pions in the background of n static baryons 
and n static antibaryons which is described by the partition function 

Z n ,n = hu -MU] n ^-mUYfeM-S[U]) exp[-/3M(n + n)]. (5.4) 

J 71/. TXi . 

The permutation factors n\ and n\ take into account the indistinguishability of 
baryons and antibaryons, respectively. Introducing the baryon chemical potential \x 
that couples to the baryon number (n — n), we obtain the grand canonical partition 
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function 

z (t i ) = J2 Z ",n e xp[(3Kn-n)] 

n,n 

= Y [VU — exp{-S[U] - (3n{M - /i) - (3n(M + /i)} 
7^ J "■ "■ 

= JvU exp{-S[U] + exp[-/3(M - //)]$[£/] + exp[-/3(M + //)]$[[/]*}. 

(5.5) 

The present calculation for baryons is consistent only for M — > oo. Hence, in order 
to obtain a non-trivial result, we also send /i to infinity keeping the difference M — /i 
finite. Then the partition function reduces to 

Z(fj) = JVU exp{-S[U] + exp[-/3(M - //)]$[£/]}. (5.6) 

In the Nf = 2 case, the nucleons transform in the fundamental representation of 
SU (2) which has a real- valued trace. Consequently, the spatial integral of the cor- 
responding Polyakov loop 

$[[/] = y rf 3 x Tr F V exp[J^ dt v 4 (x,t)\, (5.7) 

takes real values only. For Nf = 3, the trace in the fundamental representation is 
in general a complex number. However, since the baryons transform in the adjoint 
representation, the appropriate Polyakov loop averaged over spatial positions reads 

$[[/] = J d 3 x Tr A V exp[j P dt v 4 (x,t)} = J d 3 x {\Tr F Vexp[j P dt v 4 (x, t)}\ 2 - l}, 

(5.8) 

which is again real. Hence, in contrast to full QCD, both for Nf = 2 and for Nf = 3 
no complex action problem arises. Using similar methods, the thermodynamics 
of gluons in the background of static quarks have been investigated [52-54]. In 
that case, a complex action problem does arise but, at least in the Potts model 
approximation to QCD, could be solved with a meron-cluster algorithm [55]. 

It is straightforward to formulate the static baryon model on the lattice. For 
Nf = 2, the spatial integral of the Polyakov loop for a static nucleon is given by 

*[£/]= £ Tr F n^, M , (5.9) 

x t 

and, for Nf = 3, the corresponding quantity for a static baryon takes the form 

m = J2Tr A [l[V SM ] = £ {|Tr F nW " l}- (5-10) 

St x t 

Then, the lattice partition function corresponding to (5.6) reads 

Z(fi) = J] / dU x exp{-S[U] + exp[-/3(M - n)]$[U]} (5.11) 

x J SU(N f ) 
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with the pure pion lattice action S[U] given in (4.2). It should be noted that we have 
ignored the Pauli principle for static baryons occupying the same lattice site. The 
effect of this approximation is a lattice artifact that does not affect the continuum 
limit. 

The static baryon model can be used in exploratory studies of the QCD phase 
diagram at non-zero chemical potential. At /j, — 0, the pure pion lattice model has a 
finite temperature chiral phase transition, above which chiral symmetry is restored. 
For two massless flavors, this transition is second order and it is in the universality 
class of the 3-dimensional 0(4)-model. This transition extends into the (/i, T)-plane 
and is expected to turn into a first order phase transition at a tricritical point. When 
non-zero current quark masses are included, the second order phase transition line 
is washed out to a crossover and the tricritical point turns into a critical endpoint of 
the remaining line of first order phase transitions. Locating the critical endpoint in 
the (/x, T)-plane is an important issue in heavy ion collision experiments [56] which 
has been addressed using lattice methods [57]. Unfortunately, in these studies one 
is limited to small lattices due to a severe complex action problem. The small-/! 
region has also been investigated in [58,59]. One should not expect to be able to 
extract the precise location of the critical endpoint using the static baryon model 
presented here. However, as long as the model has the correct symmetry properties, 
one can at least investigate universal features of this point, which is expected to be 
in the universality class of the 3-dimensional Ising model. Moreover, one can gain 
experience with numerical methods for locating the critical endpoint in a model that 
is similar to full QCD but much simpler to simulate. 

We would like to mention that the static baryon model shares all global sym- 
metries with QCD, but is, in addition, also invariant under the intrinsic parity 
transformation P . In the Nf = 2 case, P is just G-parity and hence a symmetry of 
QCD. In the Nf = 3 case, P is not a symmetry of QCD. At the level of the effective 
theory, this symmetry is explicitly broken by the Wess-Zumino-Witten term in the 
Goldstone boson sector and by the qa term in the baryon sector. Since we have 
not included the Wess-Zumino-Witten term or the qa term in the lattice action, as 
it stands, the Nf = 3 static baryon model has the undesirable P symmetry. In a 
study of universal properties of the Nf = 3 model one may hence want to break P 
explicitly by hand in order to exactly reproduce all symmetries of QCD. 

With light up and down quarks and a very heavy strange quark, the situation 
is essentially like in the Nf = 2 case. As one lowers the strange quark mass, the 
critical endpoint is expected to move toward smaller values of /i until it reaches the 
temperature axis. For even smaller strange quark masses, the chiral phase transition 
at /j, — is first order. Lattice calculations suggest (but do not show unambiguously) 
that, for Nature's quark mass values, at /i — there is merely a crossover. In that 
case, the critical endpoint would indeed be located in the (//, T)-plane at /i > 0. 
While one cannot determine the precise location of the critical endpoint in QCD from 
the static baryon model, it should be possible to obtain qualitative (and perhaps 
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even semi-quantitative) information concerning the location of this endpoint as a 
function of the strange quark mass. 



6 Constituent Quarks on the Lattice 

In this section we use a nonlinearly realized lattice chiral symmetry in order to ad- 
dress non-perturbative questions concerning the confinement of constituent quarks. 
Again, we do not rely on the validity of a systematic low-energy expansion. The 
chiral quark model of Georgi and Manohar [36] offers an explanation for why the non- 
relativistic quark model works. However, there are some potential non-perturbative 
problems which we suggest to investigate on the lattice. 

The chiral quark model is formulated in terms of gluons, pions, and constituent 
quarks ^f(x) and ^(x) with the Euclidean action 

s[v, u, eg = J d 4 x {Mm + * 7m (^ + v , + eg* + i^n^s* 

+ ^Tr[G^]}. (6.1) 

Here G^ix) = ig s G a ^(x)\ a is the SU(N C ) gluon field, and g s is the strong coupling 
constant between gluons and constituent quarks (transforming in the fundamental 
representations of both SU(N C ) color and SU(Nf) flavor). Under gauge transfor- 
mations g(x) G SU(N C ) the gluon field transforms as 

G,(x)> = g{x)[G,{x)+d,Ux)\ (6.2) 

As usual, the field strength takes the form 

G^{x) = d^G v {x) - d v G,{x) + [G^x), G u (x)]. (6.3) 

We limit ourselves to N c > 3. The N c = 2 case is exceptional because SU(2) is 
pseudo-real and hence the pattern of chiral symmetry breaking is different [60]. 

The chiral quark model of Georgi and Manohar is based on the assumption that 
the energy scale for chiral symmetry breaking is larger than the one for confinement. 
An effective description in terms of constituent quarks which receive their mass from 
chiral symmetry breaking before they get confined by residual low-energy gluons 
should then make sense. In fact, the phenomenological success of the non-relativistic 
quark model may suggest that this picture is indeed correct. In their work, Georgi 
and Manohar provided a framework that puts the non-relativistic quark model on a 
solid field theoretical basis. In particular, their chiral quark model holds the promise 
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to share the successes of the non-relativistic quark model. For example, making 
reasonable assumptions about the confinement of constituent quarks, it describes 
baryon magnetic moments and hyperon non-leptonic decays with about 10 percent 
accuracy. In addition, in contrast to the non-relativistic quark model, the Georgi- 
Manohar model is manifestly chirally invariant and is thus automatically consistent 
with the predictions of chiral perturbation theory. 

Georgi and Manohar also pointed out some potential problems of their model. 
For example, the constituent quarks may form a pseudo-scalar bound state in ad- 
dition to the explicit pion introduced as a fundamental field. If this bound state 
is light, one would double-count the pion. Georgi and Manohar argued that this 
should not happen. The pseudo-scalar constituent quark bound state should either 
be so heavy that it lies outside the range of the low-energy theory or it might repre- 
sent the next physical pseudo-scalar state above the pion. A similar issue arises for 
the baryons. In addition to the baryons built from constituent quarks, one can form 
baryonic Skyrmions from the fundamental pion field. Again, Georgi and Manohar 
argued that double-counting should not arise because the Skyrmions may either be 
very heavy or unstable. A potentially more severe problem is related to the con- 
finement scale. In the chiral quark model the strong coupling constant is fixed by 
adjusting the color hyperfine splitting in the baryon spectrum. The estimated value 
of g s is so small that one might expect unacceptably low-lying glueball states. For 
the same reason, there might be a low-temperature deconfinement phase transition 
in the gluon sector significantly below the finite temperature chiral phase transition. 

These potential problems are impossible to address quantitatively in the contin- 
uum formulation of the chiral quark model because they involve the non-perturbative 
dynamics of confinement. In order to be able to address these issues, we formulate 
the constituent quark model on the lattice. Due to the nonlinear realization of chiral 
symmetry, global axial transformations manifest themselves as pion field dependent 
(and thus local) transformations on the constituent quark fields. In particular, the 
constituent quarks couple to the pions through the composite flavor "gauge" field. 

Using the lattice construction of a nonlinearly realized chiral symmetry presented 
in section 3, it is now straightforward to put the Georgi- Manohar model on the 
lattice. The resulting action takes the form 

S^,*,U,G\ = M£**** + ^£(**7,^^ 

X " x,n 



X,fl 



^ £ Tr[t/]£4 +A + Ul + ,U X ] - M £ Tr[MUl + M*U X ] 

^ x,n x 
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~~ 72 X! ^[^x,^Gx+(i,uG x+0 ^Gl iV + G Xtl/ G x+ o^G x+ ^ v Gl^]. (6.4) 

iis X,fJ,,U 

Here G.^ G SU(N C ) (again taking N c > 3) denotes the standard color parallel 
transporters living on the lattice links, which transform under gauge transformations 
g x eSU(N c )as 

G'x,n = 9xG x ^gl +fl . (6.5) 

In the classical continuum limit the lattice gauge field G XjIX = exp[G At (a; + /t)] is 
related to the continuum gauge field G^(x) that lives in the su(N c ) algebra. 

For gA = 0, the Dirac operator of the chiral quark model 

D X y[U, G] = M5 X ,y + ^^V X ^G X J X + ^y ~ ^ ^ G^ _ ^ 5 X - fry ) 

+ ^ {lnls A x,^x,y + l^ A x-^x,y) 

+ 2^ x ' y ~ V x,nG x ^8 x+ ^ y - K-A^l-A^-A.y)' ( 6 - 6 ) 

is 75-Hermitean, i.e. D[U,G]^ = jsD[U, G]^. Thus, its spectrum contains complex 
conjugate pairs of eigenvalues A and A*. Only the real eigenvalues are in general 
not paired. Consequently, the fermion determinant can be negative only if there is 
an odd number of negative (real and thus unpaired) eigenvalues. In the presence of 
a sufficiently large constituent quark mass M, negative eigenvalues are very much 
suppressed. Hence, the potential sign problem due to a negative fermion determinant 
is expected to be mild. A similar situation arises when one wants to simulate an 
odd number of flavors of Wilson fermions, for example, when one wants to include 
the strange quark in a QCD simulation. In this case — not directly addressing 
the potential sign problem — algorithms have been developed to perform numerical 
simulations [39-42] . It is to be expected that similar methods can be applied to the 
action presented here. Unfortunately, for g^ ^ the Dirac operator is in general 
no longer 75-Hermitean. Hence, we can no longer argue that there is only a mild 
sign problem. Only numerical simulations can tell if a severe sign problem arises 
for realistic values of g^. If this is the case, one might want to simulate at purely 
imaginary values of g^ and attempt an analytic continuation to the real physical 
value. However, in order to address the non-perturbative questions concerning the 
confinement of constituent quarks, it may be sufficient to work at gA = 0. 

In the chiral limit (M. = 0) and for r = 1 the lattice Georgi-Manohar model 
has four bare parameters: the constituent quark mass M, the axial coupling gA, the 
pion decay constant F w , and the strong gauge coupling g s . The chiral quark model 
on the lattice can be treated similarly to the low-energy baryon effective theory. In 
particular, one need not necessarily approach a second order phase transition but 
can simulate anywhere in the chirally broken phase. For example, one can tune the 
bare constituent quark mass M and the bare pion decay constant F n such that both 
the renormalized nucleon mass M£f n and the renormalized pion decay constant i^ en 
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stay at the cut-off with their ratio equal to the experimental value. Similarly, the 
bare g& can be tuned to reproduce the observed weak decay of the neutron and g s 
can be adjusted to the nucleon-delta mass splitting. With the bare parameters fixed 
in this way, one can ask if the low-energy physics in the gluon sector is consistent 
with experiment. For example, it is interesting to investigate if there is a first 
order deconfinement phase transition at an unusually low temperature or if there 
are unacceptably light glueball states. If potential problems of this kind do not 
arise, the chiral quark model may indeed be able to explain the phenomenological 
success of the non-relativistic quark model. 

Unlike the baryon effective theory, we do not expect the chiral quark model to 
represent a systematic low-energy expansion. In particular, one should not expect 
that its low-energy gluon physics is universal and depends on M, qa, F n , and g s only. 
Still, despite the fact that other details of the model are expected to influence the 
gluon sector, it is interesting to ask if the gluon dynamics can be modeled successfully 
in this way. Since one cannot expect to obtain precise results directly relevant to 
QCD, one may ask if using the machinery of lattice field theory is at all justified. 
Clearly, this is a matter of taste. However, even if only qualitative insight into the 
success of the non-relativistic quark model can be gained, this may be interesting 
enough. In addition, a lattice study of the Georgi-Manohar model would also provide 
valuable experience with dynamical fermion algorithms with exact chiral symmetry 
in a set-up that is much simpler than lattice QCD with Ginsparg- Wilson fermions. 



7 From the Chiral Quark Model to QCD? 

It is interesting to ask if a nonlinearly realized chiral symmetry can also lead to a 
new approach to lattice QCD itself. Using a linear realization of chiral symmetry 
and an explicit auxiliary pion field, similar questions have been addressed in [37, 38]. 
In this section we investigate such questions by exploring the phase structure of the 
lattice chiral quark model as a function of its bare parameters M, F n , and g s . For 
simplicity, we put qa = and r — 1. 

In the limit of infinite constituent quark mass, M = ±oo, the pions and gluons 
decouple from each other: the model splits into an SU (N c ) Yang-Mills theory and 
an SU(Nf)i, x SU(Nf)n nonlinear cr-model for the pions. The pure gluon theory 
is expected to have no bulk phase transition in the bare gauge coupling g s and, as 
a consequence of asymptotic freedom, reaches its continuum limit as g s — * 0. For 
large values of F%, the pion field U x is ordered and chiral symmetry is spontaneously 
broken to SU(Nf) L=R , while for small values of F%, it is restored. At M — ±oo 
the critical coupling F% c is independent of g s . Consequently, there is a horizontal 
line connecting the points A and B in the phase diagram of fig. 1. In the Nf = 2 
case, the phase transition is second order, while for Nf = 3, it is expected to be first 
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order. 




M 



Figure 1: Possible form of the phase diagram of the chiral quark model in the 
(M, F 2 ,l/g 2 ) -space at g& = and r = 1. Structures behind the critical surface 
connecting the points E, F , G, and H are not shown. The ranges of the variables 
are M G [-oo, oo] and F 2 , l/g 2 s G [0, oo]. 

Next we consider the limit of infinite bare coupling F 2 — oo. The pion field 
U x = u 2 x is then frozen to a space-time independent constant U . Of course, this 
does not mean that u x is a constant. This can be illustrated as follows. If U x = U 
for all x, all u x can be diagonalized by the same unitary transformation W, i.e. 
u x = Wd x W^ . The matrix d x = ds x is a constant diagonal matrix d G U(l) Nf ~ l up 
to a diagonal matrix s x G 7.{2) N f~ l obeying s 2 , = 1. The Hermitean matrices s x 
have diagonal elements ±1 and determinant 1. Hence, we can now write 

u x = Wds x W\ (7.1) 

such that indeed u 2 x = Wd 2 W^ = U is independent of x. Let us now consider the 
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other fields in the F% = oo limit. The lattice field V x ^ takes the form 

Vx,» = + u xut +fl ] = Ws x s x+p W^ = V x ^. (7.2) 

Since V x ^ is already in SU(Nf), the GL(Nf) / SU(Nf) coset decomposition yields 
Vx,n — Vx,n- Similarly, it is straightforward to show that A Xjfl = A x = A x = 0. 
We now perform the transformation V x = s x W^ on the constituent quark field 

* x = = Ws x %, U x = %V X = V x s x Wl (7.3) 

As a consequence, 

tf x ¥ x = %s x W^Ws x % = 

VaVx^x+H = ^' x S x W^Ws x S x+ ^Ws x+ ^' x+fl = <^ +A , 

V^Vl^x = % +fi s x+ ^Ws x+(l s x W^Ws x V x = (7.4) 
so that (for F% = oo) the action (6.4) takes the form 

S^', G\ = + 1 Y,(Ki»Gx,,K + » - K + ^Gl^' x ) 

X " X,[l 

~ ~l2 E ^ T [Gx^Gx+^^Gx + u^G x ^ + Gx^Gx+u^G^^^Gx^- 
iJs X,fJ,,U 

(7.5) 

Remarkably, this is nothing but the standard Wilson fermion action. Hence, in 
the limit F% = oo, g s — > 0, and M — > (point E in fig. 1), we reach the usual 
continuum limit of lattice QCD. This point is connected to the point F at strong 
coupling (g s = oo) by a line of second order phase transitions. Fine-tuning to this 
line, a massless "pion" emerges as a quark-antiquark bound state. This "pion" is 
not the Goldstone boson that results from chiral symmetry breaking, except in the 
continuum limit. Still, in Wilson's approach to lattice QCD this particle represents 
the physical pion. Then, recovering chiral symmetry requires fine-tuning as well as 
taking the continuum limit. Of course, in the chiral quark model with a nonlinearly 
realized chiral symmetry, there is a genuine Goldstone pion. However, this particle 
decouples at F% = oo. It is tempting to try to avoid fine-tuning and to represent the 
physical pion by the Goldstone boson of the nonlinearly realized chiral symmetry. 
This would suggest to approach the QCD point E from finite values of F% out of the 
chirally broken phase. As one approaches the point E, one should then avoid the 
critical surface connecting the points E, F, G, and H because otherwise one would 
double-count the pion. However, it is not clear if one will obtain a new approach to 
QCD in this way. For example, it is likely that the explicit Goldstone pion decouples 
and that the QCD pion must still be assembled as a quark-antiquark bound state. 
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Let us now consider the g s = plane in the phase diagram. Then the gluon field 
is completely frozen to a pure gauge G X)fl = g x g x+ n- It should be noted that the 
continuum limit of lattice QCD is taken by sending g s — > in a controlled way, not 
by simply putting g s = 0. In particular, at g s = there is no gluon field and thus no 
confinement. In this case, we transform the fermion field as ty x — g x ^' x , ty x = ^ x g x) 
so that 

**** = Kd9x% = %%, 

^ x G x ^ x+ix = %gtg x gl +fl g x+i >y x+ii = 

Vx+fiGlflx = K+agl+agx+nglgxK = (7.6) 

Hence, for g s = the action (6.4) reduces to 

X ^ X,fl 



?2 



- jE nulU x+(l + Ul +fl U x ] - M £ Tr[^C/t + 

(7.7) 

which has exactly the same form as the baryon action (3.22). Of course, here the 
field ^ x describes a colored constituent quark. Since we have now switched off the 
strong interactions, these quarks are no longer confined. 

For Nf = 2, a possible form of the phase diagram in the (M, F^)-plane (at g s = 0) 
is depicted in fig. 2. It also includes the Aoki phase not shown in fig. 1. In the limit 
of infinite constituent quark mass, M = ±oo, the model reduces to an SU(2)l <g> 
SU(2)r = 0(4) nonlinear a-model for the pions. For F% > F% c , the pion field U x is 
ordered and hence chiral symmetry is spontaneously broken to SU(2)l = r = 0(3), 
while for F% < F% c , it is restored. It is interesting to consider negative values of F% 
as well (still keeping M = ±oo). On a bi-partite lattice one can map the theory back 
to the one with positive values of F% by changing the sign of U x on all odd lattice 
sites. Hence, for F% < —F% c , there is an "antiferromagnetically" ordered phase in 
which chiral symmetry is again broken from 0(4) to 0(3). Hence, for M = ±oo, 
there are two phase transitions, one at F% c separating the usual "ferromagnetic" 
chirally broken from the "paramagnetic" chirally symmetric phase and another one 
at — F% c separating the "paramagnetic" from the "antiferromagnetic" broken phase. 

At M = and F% = oo (point E in figs. 1 and 2) there is a second order phase 
transition at which the constituent quark mass vanishes. In addition, there are un- 
physical second order phase transitions at M + 2nr = (with n G {1,2,3,4}) at 
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Figure 2: Possible form of the phase diagram of the chiral quark model in the 
(M, F^) -plane at g s — 0, gA — 0, and r — 1. The ranges of the variables are 

M,F*e [-00,00]. 

which some of the doubler fermions become massless. In the context of Wilson's 
lattice QCD the region of M < was explored by Aoki [61]. He identified a rich 
structure including a phase in which parity and flavor are spontaneously broken. It 
should be noted that for g& = the model (7.7) has a symmetry under the replace- 
ment of M by — (M + 4r). This replacement can be compensated for by changing 
the sign of ^ x on all even sites and the sign of ^ x on all odd sites. This symmetry 
manifests itself in the phase diagram of fig. 2. Our sketch of the Aoki phase in 
fig. 2 is inspired by its form in Wilson's lattice QCD. The precise structure of the 
phase diagram can only be inferred from numerical simulations. The phase diagram 
of fig. 2 is a possible interpolation between the limits M = ±00 and F% = 00, 
partly inspired by the phase diagram of the Smit-Swift model and other related 
Higgs- Yukawa models [62-64]. Attempts to construct chiral gauge theories as the 
continuum limit of such lattice models have failed because it was impossible to re- 
move the doubler fermions, both at weak and at strong Yukawa couplings [65] . Such 



30 



problems do not affect our construction of lattice effective theories with nonlinearly 
realized chiral symmetry because, in that case, there is no need to approach a second 
order phase transition to extract low-energy continuum physics. However, attempts 
to construct full QCD with a manifest nonlinearly realized chiral symmetry on the 
lattice are likely to fail for the same reasons as the Smit-Swift model did. Still, it is 
interesting to ask if the lattice chiral symmetry of Ginsparg- Wilson fermions can be 
realized nonlinearly. We leave such questions for future investigations. 



8 Conclusions 

In this paper we have constructed lattice theories with baryon or constituent quark 
fields on which chiral symmetry is nonlinearly realized. In this framework the 
fermion doublers are removed by the Wilson term without breaking chiral symmetry 
explicitly. This is possible because fermion fields do not contribute to anomalies in 
theories with nonlinearly realized chiral symmetries. Instead the anomalies must 
be included explicitly in the form of Wess-Zumino-Witten and Goldstone-Wilczek 
terms. The lattice construction of such terms remains an interesting challenge for 
future investigations. 

The use of nonlinearly realized chiral symmetry is well established in the context 
of chiral perturbation theory which provides a power counting scheme for a system- 
atic low-energy expansion. However, also non-perturbative problems may arise in 
the low-energy regime. Some of these problems can be dealt with analytically [28]. 
In general, however, such problems require some numerical work. While one can 
still proceed without Monte Carlo methods in investigations of few-nucleon systems 
[30,31], this is likely to change as one moves on to larger nuclei and ultimately to 
nuclear matter [32]. Our lattice formulation provides a framework in which future 
numerical studies of this kind can be performed on a solid theoretical basis. 

Even if one leaves the framework of a systematic low-energy expansion, the 
nonlinear realization of chiral symmetry on the lattice may provide insight into 
interesting physical questions. For example, as we have seen, in a model of static 
baryons at non-zero chemical potential one can investigate universal aspects of the 
QCD phase diagram, like the nature of the critical endpoint of the chiral phase 
transition line. In the static baryon approximation, one cannot expect to obtain 
results directly relevant to QCD for non-universal quantities such as the location 
of the critical endpoint. Still, one may be able to collect, for example, qualitative 
information on how the critical endpoint moves as the strange quark mass is varied. 

Another interesting example is the chiral quark model of Georgi and Manohar. 
It holds the promise to explain the success of the non-relativistic quark model and 
has successfully described a variety of physical properties of hadrons. However, 
there is a number of potential dynamical problems that can only be addressed non- 
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perturbatively. Our lattice formulation allows us to investigate if these potential 
dynamical problems do arise or not. In particular, the confinement of constituent 
quarks — out of reach of conventional methods — can be directly addressed with lat- 
tice methods. These investigations will allow one to decide if the chiral quark model 
indeed provides a satisfactory explanation for the success of the non-relativistic 
quark model. Along the way, one will be able to sharpen the algorithmic tools of 
dynamical fermion simulations in a model that shares many features with QCD - 
most important a manifest chiral symmetry — but which should be much simpler 
to simulate than Ginsparg- Wilson fermions in the chiral limit. 

At this stage, it is not clear if lattice theories with a nonlinearly realized chiral 
symmetry can also provide a formulation of full QCD beyond low-energy effective 
theories. In order to answer this question, one may attempt to relate lattice QCD 
with Ginsparg- Wilson fermions to a lattice theory with explicit pion fields and a 
nonlinearly realized chiral symmetry. If such a connection can be found analyti- 
cally, one would expect that the anomalies of Ginsparg- Wilson fermions manifest 
themselves as Wess-Zumino-Witten and Goldstone-Wilczek terms of the lattice pion 
field. This is an interesting project for future studies which holds the promise to 
further deepen our understanding of the non-perturbative regularization of chiral 
symmetry. 
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